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(57) Abstract 

A method for process- 
ing an image from a multicore 
fibre or multifibre, compris- 
ing the steps of (1) calculat- 
ing, for each core of the mul- 
ticore fibre or the multifibre, 
an intensity equal to the sum 
of the intensities of the image 
pixels located inside a prede- 
fined area defining said core; 
(2) calculating a standardised 
intensity for each core by di- 
viding the transmitted inten- 
sity thereof, calculated as in 
step 2, by the intensity trans- 
mitted by the same core in a 
reference image; and (3) cal- 
culating an interpolated value 
of the 'intensity of each image 
pixel depending on the stan- 
dardised intensities of three cores in a predetermined triangular mesh, and on the position of the pixel relative to said three cores. The 
method is useful in medical analyses. 
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(57) Abr^ 

L' invention concerne un proc6d6 de traitement d'une image obtenue par une fibre multicoeur ou une multifibre, comportant les 6tapes 
suivantes: 1) pour chaque coeur de la fibre multicoeur ou de la multifibre, calcul d'une intensity, 6gale a la somme des in ten sites des pixels 
de 1' image qui sont situ6s a rinterieur d'une zone, pr6d61imit£e, de delimitation dudit coeur; 2) calcul d'unc intensity normalises de chaque 
coeur, en divisant son intensity transmise, calculee comme dans l'ltape (2), par V intensity transmise par ce mcmc coeur dans une image de 
reference; 3) calcul d*une valeur interpolee de Pintensite* de chaque pixel de I 'image, en fonction des intensitts normalis^es de trois coeurs 
d'un maillage triangulaire pr6d6termine\ et de la position relative du pixel par rapport a ces trois coeurs. Application a V analyse meclicale. 
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1 

PROCEDE DE TRAITEMENT D ' IMAGES OBTENUES PAR FIBRES 
MULTICOEURS OU MULTI FIBRES , EN PARTI CULIER D » IMAGES 

ENDOSCOPIQUES 

DESCRIPTION 

Domaine technique et art anterieur 

L 1 invention concerne le domaine du traitement 
d' images, notamment le domaine de traitement d ' images 
obtenues en microendoscopie , par fibres optiques 
multicoeurs . 

Une fibre multimode classique comporte, comme 
illustre sur la figure 1A, un coeur 1 et un manteau 3. 

Une fibre multicoeurs est un faisceau de 
fibres, fondu et etire, qui forme done un ensemble 
continu. Le manteau de chaque fibre individuelle st 
fondu avec les manteaux des coeurs voisins . 

A l'interieur d 1 une fibre multicoeurs on ne 
peut distinguer que des coeurs individuels, le manteau 
des fibres etant devenu en quelque sorte collectif. 

La figure IB represente une vue en coupe 
transversale d ' une fibre multicoeurs, les coeurs 24 et 
les manteaux 26 etant groupes a l'interieur d ' une 
premiere gaine 28, par exemple en silice, et d • une 
25 seconde gaine 30, dite gaine externe ou revetement 
"noir". Le diametre exterieur D : de 1' ensemble peut 
etre par exemple de l'ordre de 200 a 500 um. 

La figure 1C est une vue agrandie de la portion 
32 du faisceau de coeurs. Sur cette figure 1C, il 
apparait que les coeurs ont des sections transversa les 
de forme variable, plus ou morns homogene . En 
particulier, le diametre d de chaque coeur, e'est-a- 
dire la plus grande distance separant deux points d 1 un 
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meme coeur, varie d'un coeur a 1* autre. Typiquement d 
peut, par exemple, varier entre 3 et 4 vim pour une meme 
fibre multicoeurs. De meme, la distance moyenne d'un 
coeur a 1* autre n'est pas uniforme et peut par exemple 

5 varier, pour une meme fibre multicoeurs, de 3 a 3,5 prn. 

La notion de fibre multicoeurs est a distinguer 
de celle de multifibres, qui est un assemblage ou un 
faisceau de fibres independantes placees con jointement 
et 6ventuellement collees en bout. L' invention 

10 s* applique egalement aux multifibres. 

Les fibres multicoeurs et les multifibres sont 
utilisees en imagerie, notamment dans le domaine 
medical. L ' endoscopie, et en particulier la 
microendoscopie, permet au praticien d'acquerir des 

15 informations, ou images, des parties interieures au 
corps humain, telles que l'estomac, les poumons ou le 
coeur . 

Un dispositif pour la mise en oeuvre d'une 
telle technique est represents schemat iquement sur la 

20 figure 2, ou la reference 2 designe une source de 
lumiere qui est focalisee par une lentille 4 a 1' entree 
d'un guide de lumiere 6. Ce dernier est en fait le plus 
souvent relie a une pluralite de fibres optiques 8, 10 
disposees a la peripheric d 1 une fibre multicoeurs 12. 

25 Un faisceau d'eclairage 14 peut ainsi etre dirige sur 
une zone 16 d'un objet ou d J un organe a observer, qui 
reflechit un rayonnement 18 a 1' entree 20 d'une fibre 
multicoeurs 12. Cette derniere comportant un faisceau 
coherent de coeurs individuels, ceux-ci transmettent 

30 done la lumiere de maniere ordonnee entre eux, et 
1 ' image obtenue en sortie 22 de la fibre multicoeurs 
correspond a 1 1 image formee a 1' entree 20. Des moyens 
pour memoriser, analyser et/ou representer I ' image 
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peuvent etre egalement prevus en combinaison avec ce 
dispositif . 

Cette technique d'imagerie est decrite par 
exemple dans les articles de A. Katzir : "Optimal 
Fibers in Medicine", Scientific American, vol.260 (5), 
p. 120-125, 1989 et "Optimal Fiber Techniques 
(Medicine)", Encyclopedia of Physical Science and 
Technology, vol. 9, p. 630-646, 1987. 

Pratiquement , une fibre multicoeurs telle que 
la fibre 12 peut comporter environ 700 a 10000 coeurs, 
pour des applications a la microendoscopie . 

Dans l 1 image obtenue par un microendoscope, on 
constate la presence d 1 un reseau des coeurs du faisceau 
de fibres. Ce reseau est constitue de regions presqjje 
circulaires, dont chacune correspond a un coeur. 

Le praticien ou la personne qui examine une 
image obtenue par microendoscopie, ou, plus 
generalement, obtenue a l'aide d'une fibre multicoeurs, 
ne peut done utiliser 1 ' image brute, directement 
obtenue en sortie de la fibre. Un traitement numerique 
de cette image est necessaire. 

L' article de R. Conde et al . Intitule 
"Comparative Measurements of Image Quality in Image 
Guides" paru dans Biomedical Optics' 94 decrit un 
traitement d'une image test d'une mire de Ronchi 
obtenue par fibre multicoeur. Ce document ne donne 
cependant aucune precision sur la mise en oeuvre 
pratique, qui pose de gros problemes en particulier du 
fait de la masse importante de donnees a traiter, 
presentes dans chaque image obtenue par une fibre 
multicoeur (typiquement de l'ordre de 512x512). 

Expose de l'invention 
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L ' invention a pour objet un procede de 
traitement d'une image digitalisee obtenue par une 
fibre multicoeur, ou par une multifibre, comportant : 

1. pour chaque coeur de la fibre multicoeur ou de la 
5 multifibre, calcul d'une intensite, egale a la somme 

des intensites des pixels de 1 1 image qui sont situes 
a l'interieur d'une zone, predel imitee , de 
delimitation dudit coeur 

2. calcul d'une intensite normalisee de chaque coeur, en 
10 divisant son intensite transmise, calculee comme dans 

l'etape 1, par 1 1 intensite transmise par ce meme 
coeur dans une image de reference 

3. calcul d'une valeur interpolee de l 1 intensite de 
chaque pixel de 1' image, en fonction des intensites 

15 normalisees de trois coeurs d'un maillage 

triangulaire predetermine et de la position relative 
du pixel par rapport a ces trois coeurs. 

Ce traitement permet, dans 1' image obtenue, 
d'eliminer la presence du reseau des coeurs du faisceau 
20 de fibres de la fibre multicoeur ou de la multifibre. 
Un tel reseau est constitue de regions presque 
circulaires, dont chacune correspond a un coeur. Ainsi, 
la distribution d f intensite de 1' image devient 
continue. 

25 L'etape de normalisation, ou de calibrage de la 

reponse de chaque coeur, permet de tenir compte des 
differences de reponse des coeurs a un meme signal. Ces 
differences de reponse sont dues entre autres au fait 
que la distance entre les coeurs varie dans le reseau 

30 de fibres, ce qui a pour consequence des couplages 
entre les coeurs : pour chaque mode de propagation a 
l'interieur d'un coeur individuel, il existe un champ 
evanescent reparti spatialement dans les zones 
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intercoeurs et dans les coeurs voisins. Ce champ est 
d'autant plus important que l'ordre du mode est eleve. 

II en resulte que, pour une intensite incidente Iq a 

l f entree d'un coeur, celui-ci transmet en sortie une 

5 certaine intensite Iq-Iq, tandis que les coeurs voisins 

transmettent une intensite ±q. Des variations de 

.distance entre les coeurs voisins entrainent done des 

variations d'intensite ainsi transmises i 0 d'un coeur 

voisin a un autre coeur voisin. 
10 Enfin, la methode d ' interpolation a partir d'un 

maillage triangulaire predetermine est particulierement 
bien adapte pour obtenir une image continue, dans vw le 
cas d'une image obtenue par une fibre multicoeur ou par 
une multifibre. 

15 La zone predeterminee autour de chaque coeur 

peut etre definie sur 1 1 image de reference, cette 
derniere etant transposee sur 1 ' image objet apres .que 
celle-ci ait ete obtenue. Ainsi, il suffit de faire une 
seule fois la determination des zones de delimitation 
20 autour de chaque coeur.' On utilise ensuite, pour chaque 
image, la meme zone de delimitation. 

La zone de delimitation autour de chaque coeur 
peut etre la surface delimitee, dans 1 1 image de 
reference, par les mediatrices passant : 
25 - entre le barycentre dudit coeur et le barycentre de 
chacun des coeurs voisins, 
- ou, entre le pixel d'intensite maximum de chaque 
coeur et le pixel d'intensite maximum de chacun des 
coeurs voisins . 
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La premiere definition (a 1 * aide des 
barycentres de coeur) est tres pratique car elle ne met 
en oeuvre que des considerations geometriques . 

En ce qui concerne l'etape d ' interpolation de 
5 la valeur de l'intensite de chaque pixel de 1 ' image 
objet, chaque pixel peut etre avantageusement 
prealablement muni de 1 • identification des trois coeurs 
du maillage triangulaire qui lui correspondent. La 
encore, cette identification peut etre faite des le 
10 depart, sur 1 ' image de reference, et etre ainsi 
disponible par la suite pour toute image objet. 

L' image objet interpolee peut ensuite etre 
rehaussee. Cette etape permet d'ameliorer le contraste 
de 1 ' image . 
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Le procede tel que decrit ci-dessus peut 
egalement comporter les etapes prealables suivantes de 
traitement d 1 une image digitalisee de reference : 
1'. recherche du barycentre, ou du pixel d'intensite 

maximale, de chaque coeur dans 1 ' image de reference. 
2'. calcul de la zone de delimitation autour de chaque 

coeur . 

Ces etapes permettent de determiner, une fois 
pour toutes, les elements de reference (point a 
l'interieur de chaque coeur : barycentre ou pixel 
d'intensite maximale ; zone de delimitation autour de 
chaque coeur, qui pourra etre reutilisee par la suite 
dans tout traitement d' image objet). 

Avantageusement, on realise, apres 

digitalisation des pixels de 1 ' image de reference, le 
calcul d'un masque numerique qui delimite la zone de 
1' image a l'interieur de laquelle, pour chaque coeur, 
la recherche du barycentre, ou du pixel d'intensite 
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maximale, et le calcul de la zone de delimitation 
doivent etre effectues. Cette etape de calcul d'un 
mas.que numerique permet, par la suite, d'appliquer les 
operations de traitement uniquement dans la zone 
delimitee par ce masque : on limite ainsi le nombre 
d' operations et le temps de calcul. 

La recherche du barycentre, ou du pixel 
d'intensite maximale, peut etre faite apres avoir 
determine, pour chaque pixel de 1 • imaqe de reference, 
1' ensemble des pixels qui 1 ' entourent et qui sont 
d'intensite superieure a un certain seuil donne T . On 
procede ainsi par reqroupement , a partir de chaque 
pixel, et on attribue un numero au coeur correspondant , 
c'est-a-dire a 1 ' ensemble des pixels ainsi regroup.es 
pour chaque coeur. 

Pour chaque pixel de 1 1 image de reference, on 
peut determiner, notamment a l'aide de ses coordonnees, 
si on est en presence d'un pixel situe au bord de 
1' image, ou pas. 

20 Ainsi, les determinations des pixels 

d'intensite superieure a T peut evoluer lineairement , 
et changer de direction des qu'un pixel d'intensite 
inferieure a un seuil T est rencontre, ou bien des 
qu'un pixel de bord est detecte. 
25 Dans 1' image de reference, la zone de 

delimitation associee a chaque coeur peut etre obtenue 
en associant, au barycentre ou au pixel d'intensite 
maximale de chaque coeur, une surface limitee par les 
mediatrices entre ledit coeur et chacun des coeurs 
voisins. Ainsi, on attribue chaque pixel a une zone 
autour d'un coeur. En outre, le calcul d'une mediatrice 
entre deux points de coordonnees definies est aise a 
realiser. Le calcul des surfaces delimitees par les 
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mediatrices est done lui aussi aise a realiser, Cet 
aspect est important, dans la mesure ou le nombre de 
coeurs a traiter, & l'interieur d 1 une meme fibre 
multicoeur, peut etre tres eleve (de 1 1 ordre de 
5 plusieurs milliers) - 

Dans le cas ou les coeurs fusionnes sont 
presents, ce qui peut arriver tres souvent, mais avec 
une importance variable d 1 une fibre multicoeur a 
1' autre, une etape de determination de la presence de 
10 coeurs fusionnes et du barycentre ou du pixel 
d'intensite maximale, de chacun des coeurs de chaque 
groupe de coeurs fusionnes peut etre realisee. Ainsi, 
tous les coeurs, sans exception, y compris ceux qui 
sont fusionnes avec des coeurs voisins, sont traites de 
15 la meme mani£re dans 1' image de reference, et done 
egalement dans toute image objet. 

Avantageusement, la sous etape de traitement 
d 1 un groupe de coeurs fusionnes comporte : 

la determination du barycentre geomet rique des pixels 
20 faisant partie du groupe de coeurs fusionnes. 

la determination d'une direction de fusion, d'une 
direction de coupure, et de points de coupure du 
groupe de cours fusionnes. 
- la coupure du groupe de coeurs fusionnes selon la 
25 direction de coupure, en chaque point de coupure. 

Avantageusement, le procede, avec etape 
prealable de traitement de 1 1 image de reference, peut 
comporrer en outre une etape prealable de determination 
d ! un maillage triangulaire de 1 ' image de reference, 
30 telle que chaque pixel fasse partie d'un triangle. 

Ainsi, pour toute image objet, on dispose d'un 
maillage triangulaire predetermine permettant, pour 
tout pixel de cette image objet, de calculer une 
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interpolation en fonction des intensites normalisees 
des pixels constituant les sommets du triangle dont 
fait partie le pixel considere de 1 ' image objet. 

L'etape de determination d'un maillage 
triangulaire peut comporter : 

- une sous-etape de determination, pour chaque 
barycentre ou point d'intensite maximale, des 
barycentres ou points d' intensites maximales 
immediatement voisins, et de determination des pixels 
qui appartiennent a chaque triangle forme par deux 
barycentres ou points d'intensite maximale, choisis 
voisins l'un de 1 ' autre parmi les barycentres ou 
points d'intensite maximale immediatement voisins du 
barycentre ou point d'intensite maximale consider^, 
et par ce dernier 

- pour les pixels n'ayant aucun triangle attribue par 
la sous-etape precedente, une sous etape de : 

* recherche des deux pixels les plus proches 
appartenant a un triangle, 

* determination des trois barycentres, ou points 
d'intensite maximale, les plus proches du pixel 
considere, parmi les six barycentres ou points 
d'intensite maximale, constituant les sommets des 
deux triangles auxquels appartiennent les deux 
pixels les plus proches du pixel considere. 

Le fait de distinguer entre les triangles ayant 
deux sommets alignes horizontalement et les autres 
permet de reduire le temps de calcul par un facteur 
environ egal a 2. 

La technique de traitement d' images objet de 
1 ' invention peut etre appliquee a des images en noir et 
blanc ou bien a des images couleurs. Dans ce dernier 
cas, 1' image est decomposee en trois images de base, 
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une pour chacune des couleurs Rouge, vert, Bleu (R, G, 
B) . Chaque image de base peut alors etre traitee par le 
procede objet de la presente invention. Pour fournir 
une image en couleurs, le dispositif correspondant fait 
5 usage d'une camera adapt6e. 



Breve description des figures 

De toute fagon, les caracteristiques et 
avantages de 1' invention apparaitront mieux a la 
10 lumiere de la description qui va suivre. Cette 
description porte sur les exemples de realisation, 
donnes a titre explicatif et non limitatif, en se 
referant a des dessins annexes sur lesquels : 

- les figures 1A, IB et 1C representent des 
15 vues en coupe transversale d'une fibre optique et d'une 

fibre mult icoeur , 

- la figure 2 est une representation 
schematique d'un dispositif d ' endoscopie , 

- la figure 3 est une representation 
20 schematique, en bloc, du systeme informatique utilise, 

- les figures 4A et 4B representent des images 
brutes obtenues a l'aide d'une fibre multicoeur, en 
lumiere blanche (image de reference : figure 4A) et 
pour une image objet (figure 4B) , 

25 - la figure 5 represente des coeurs d'une image 

obtenue, des pixels de l 1 image, et des zones 
geometriques de 1 1 image permettant de traitement de 
celle-ci, 

- la figure 6 illustre schemat iquement un 
30 procede d : interpolation d'une image obtenue, 

- les figures 7A a 7E illustrent des etapes de 
traitement d'une image de reference, 
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- la figure 8 represente une partie d'une 
image, avec des coeurs fondus, 

- les figures 9A a 9C representent des etapes 
de traitement des groupes de coeurs fusionnes, 

5 - la figure 10 est un organigramme resumant des 

etapes du traitement d'une image de reference 
(procedure de segmentation) , 

- les figures 11A a 11C schematisent des etapes 
de determination d 1 un maillage t riangulaire , 

10 - la figure 12 illustre une methode d'extension 

d 1 histogramme , pour realiser une etape de rehaussement 
de 1 ' image . 

Description detaillee de modes de realisation -ide 

15 1 ' invention 

Le traitement d 1 image endoscopique , conforme- a 
1* invention, peut etre mis en oeuvre a i'aide d'un 
dispositif tel que celui illustre sur la figure 2. 
Cette figure a deja ete partiellement commentee ci- 

20 dessus. Les images de l'objet 16 obtenues en sortie 22 
peuvent etre visualisees par une camera ou tout autre 
dispositif de visualisation 34, par exemple relie a un 
ordinateur 36 comportant un ecran de visualisation 37 
et un clavier 38. L 1 ecran permet en particulier de 

25 representer une image apres un traitement conforme a 
1 1 invention . 

L'ordinateur comporte une section de. calcul 
avec microprocesseur et toutes les composantes 
electroniques necessaires au traitement des images 
30 obtenues. 

La figure 3 est une representation simplifiee, 
en bloc, des composantes de l'ordinateur. Un 
microprocesseur 39 est relie, par un bus 40, a des 
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memoires RAM 41 pour stocker cies donnees, a une memoire 
ROM 42 dans laquelle les instructions du programme de 
traitement realise peuvent etre memonsees. Le systeme 
comporte en outre le dispositif de visualisation, ou 
5 ecran, 37 et le clavier 38 . 

Un ordinateur personnel, du type POWER PC (800- 
80 AV) peut etre utilise en tant qu 1 ordinateur 36 dans 
le cadre de la presente invention. 

Le schema represents sur la figure 2 permet une 

10 exploitation " immediate" , par un operateur, des images 
visualisees. Dans le cas d' images medicaies, ledit 
operateur peut etre un praticien qualifie, pouvant 
analyser des images traitees. 

Dans d'autres cas, 1* operateur, par exemple un 

15 praticien generaliste, ne saura pas analyser 1 1 image 
traitee. Dans ce cas, cette derniere peut etre 
memorisee et envoyee a un site d 1 analyse, par exemple a 
un laboratoire d 1 analyse specialise. Ce dernier 
etablit, au vu de 1' image, un diagnostic, qu'il 

20 communique ensuite au praticien generaliste. Celui-ci 
peut alors par exemple diriger, en fonction du 
diagnostic, le patient vers un medecin specialise. 
L' ensemble du dispositif illustre sur la figure 2 
(fibre multicoeur, camera 34, ordinateur 36) presente 

25 un encombrement faible, et peut done facilement faire 
partie de 1 1 instrumentation dont peut disposer un 
medecin generaliste. La transmission des images, a 
distance, au site d 1 analyse ou au laboratoire d'analyse 
peut etre faite par tout moyen approprie, par exemple 

30 en stockant les images sur des disquettes envoyees 
audit laboratoire. Un autre moyen consiste a envoyer 
les donnees au site ou au laboratoire d'analyse par un 
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reseau electronique de type Internet auquel le 
generaliste et le laboratoire sont relies. 

Les images visualisees sur l'ecran 37 sont les 
images obtenues en sortie 22 de fibres multicoeurs et 
sont du type de celles representees sur les figures 4A 
et 4B. Chaque image represente un reseau des coeurs 
individuels de la fibre multicoeur ou de la multifibre. 
Ce reseau est constitue de regions presque circulaires, 
dont chacune correspond a un coeur. Sous sa forme 
brute, 1' image ainsi directement obtenue est 
inutilisable par un praticien ou un operateur, et il 
est necessaire de la traiter pour que celui-ci puisse 
la lire et 1 1 interpreter . 

Dans tous les cas, le procede selon 1' invention 
permet au praticien ou a 1' operateur de retrouver 
pratiquement, a partir de 1' image non traitee et 
inutilisable, une image pouvant etre utilisee. 

La figure 4A represente schemat iquement 1 1 image 
obtenue en exposant la fibre multicoeur a un eclairage 
blanc uniforme : Pour cela, une source de lumiere 
blanche est disposee a 1' entree de la fibre multicoeur, 
a la place de l'objet 16 sur la figure 2. 

La figure 4B represente 1 ■ image obtenue 
lorsqu'un corps 16 (figure 2) est examine : la source 2 
emet un rayonnement qui va eclairer ce corps, qui lui- 
meme renvoie un rayonnement dans la fibre multicoeur ; 
ce rayonnement reflechi constitue 1 • image comportant 
1 ' information qui va etre traitee. 

Chacune des images (image de reference, image 
objet) est en fait constitute d'un ensemble de NxN 
points, ou pixels. A chaque pixel correspond une valeur 
numerique, I'intensite digitalisee, qui est une 
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fonction de la position p(x,y) du point, ou pixel, dans 
le plan de 1 1 image . La figure 5 represente la 
superposition d'une partie de 1 1 image brute, 
representant les images des coeurs, avec des pixels de 
5 1' image, symbolises par des croix sur le reseau de 
1' image des coeurs. Sur cette figure 5, on voit que les 
pixels sont situes a 1 ' interieur des coeurs, d'autres 
pixels etant situes entre des coeurs voisins. Le nombre 
de pixels et le nombre de coeurs representee sur la 

!0 figure 5 ne sont pas dans un rapport qui correspond a 
la realite. En pratique, 1 1 image peut comporter la 
trace de plusieurs milliers de coeurs, pour plus de 
250000 pixels (512x512) . 

Par consequent, la partie de 1 ' image qui 

15 correspond a chaque coeur est formee d'un ensemble de 
pixels ayant certaines intensites. II est possible de 
schematiser chaque coeur par un pixel unique qui 
rassemble 1 ' information, ou intensite, moyenne des 
pixels qui forment son image. Ainsi, comme on le verra 

20 plus loin, on peut associer a chaque coeur les 
coordonnees d'un point qui est par exemple le 
barycentre des pixels rassembles dans son image. Ce 
barycentre est affecte d'une intensite pouvant etre la 
somme des intensites des pixels qui forment 1 ' image du 

25 coeur considere. 

Pour eliminer 1 1 image du reseau des coeurs, il 
faut determiner, pour les pixels situes dans les 
regions intercoeurs, de nouvelles valeurs d ! intensite. 
Celles-ci sont obtenues par interpolation des 

30 intensites des coeurs voisins. Selon la presente 
invention, un traitement a partir des trois coeurs "les 
plus proches" de chaque pixel permet une interpolation 
efficace dans un temps acceptable du point de vue 
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pratique. Pour chaque pixel, les trois coeurs 
permettant 1 ' interpolation, ou le triangle delimite par 
ces. trois coeurs, peuvent etre reperes sur 1 * image de 
reference : dans la mesure ou la meme fibre multicoeur 
est conservee d ' une image objet a une autre, la meme 
image de reference peut etre etablie, une fois pour 
toutes, pour toutes les images objets. 

Ainsi, si la figure 5 represente une partie de 
1' image de reference, des triangles d ' interpolation 54, 
56, 58 seront identifies. Chaque pixel, a l'interieur 
d'un triangle donne, se verra affecte d ' un code 
representant ledit triangle. Un tel code peut etre par 
exemple les numeros des coeurs, si ceux-ci sont 
numerotes, ou bien les coordonnees des centres ou des 
barycentres de ces coeurs. La premiere solution est 
preferable car elle necessite beaucoup moins de place 
pour' stocker 1 ' information relative an codage. 

II a ete dit plus haut que chaque coeur pouvait 
etre repere par un seul point, affecte d ' une intensite 
globale. En fait, pour avoir une evaluation exacte de 
la reponse de chaque coeur, on associe a chaque coeur, 
conformement a la presente invention, la somme des 
intensites des pixels de 1 ' image qui sont situes a 
l'interieur d • une zone 52, dite zone de delimitation 
dudit coeur. Une technique de delimitation de ces zones 
sera expliquee plus loin. 



La procedure d ' interpolation peut etre 
expliquee en liaison avec la figure 6, qui represente 
30 1' intensite I (coordonnee 2) integree (c ' est-a-dire 
obtenue comme il vient d'etre explique) associee a 
chaque coeur, en fonction de la position (XY) du coeur, 
ou de son barycentre, dans 1 1 image . Sur la figure 6, on 
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s'interessera plus particulierement aux coeurs reperes 
par les lettres a, b, c dont les coordonnees sont 

respectivement (X a , Y a ), (X b , Y b ) , <X C , Y c ) et aux 

points A(X a/ Y a , I a ) , B(X b/ Y b/ I b ) , C(X C/ Y c , I c ) cians 

5 l'espace XYZ. L 1 interpolation bilineaire est calculee 
sur la surface triangulaire determinee par les 
barycentres des trois coeurs : les intensites des 
pixels p contenus dans la zone triangulaire entre les 
coeurs a, b, c (zone 60 sur la figure 6) sont 

10 remplacees par des valeurs d'intensite interpolee. La 
valeur d'intensite interpolee, pour chaque pixel p, est 
fonction de la distance de ce pixel aux sommets a, b, c 
ainsi que des intensites affectees a chacun de ces 
derniers . 

15 De maniere plus precise, si on considere le 

point interpole P dans le triangle A, B, C, on a : 

0P=0A + AP=0A + (aAB + PAc) , 

ou, en terme de coordonnees cartesiennes : 

X p =X a +a ( X b -X a ) +p ( X c -X a ) (2.1) 

20 Y p =Y a +a ( Y b - Y a ) + P ( Y c »Y a ) (2.2) 

Z p =Z a +a ( Z b - Z a ) +P ( Z c - Z a ) (2.3) 

De ces equations, on peut deduire Z p , qui est 

la valeur d'intensite interpolee associee au point p. 

Par rapport aux techniques de traitement 

25 d' images par convolution, la presente methode ne 
necessite que de realiser que deux multiplications et 
trois additions par point ou pixel. Une technique de 
convolution necessite, pour chaque pixel, au moins N 
multiplications et N additions, avec N>2 ou 3 (le plus 

30 souvent N>3 ) . 
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Un traitement pouvant etre applique a 1 • image 
de reference va maintenant etre decrit . Ce traitement 
permet, pour chaque coeur individuel, d' identifier un 
5 point particulier dans ce coeur, de lui affecter une 
intensite cor respondant a 1' intensite globale vehiculee 
par le coeur (cette partie du traitement est resume sur 
la figure 10 sous forme d 1 organigramme ) , et, pour 
chaque pixel, de reperer les trois coeurs, qui 

10 permettront de calculer son intensite interpolee de la 
maniere decrite ci-dessus. 

Tout d'abord, il est preferable de delimiter la 
surface circulaire du reseau des coeurs qui contient 
1 ' information utile pour le traitement ulterieur. Cette 

15 delimitation permet un gain du temps de calcul et du 
nombre d 1 operations . 

Prat iquement , un masque numerique pour 1 ' image 
de la fibre est defini, qui permet de distinguer d ' une 
part les pixels de la fibre multicoeurs, et d' autre 

20 part les pixels situes dans la region entre la fibre 
multicoeurs et la gaine externe de la fibre (gaine 
referencee 30 sur la figure IB) et au-dela de cette 
region. La figure 7A represente 1 1 image de reference 
obtenue en eclairage blanc. La reference 62 designe la 

25 surface circulaire du reseau des coeurs de la fibre 
multicoeurs : cette surface contient 1 1 information 
utile pour le traitement de l' image. La determination 
d'un masque numerique , qui delimite cette surface, 
permet d'appliquer les operations ulterieures seulement 

30 dans cette region, et done de reduire le nombre 
d'operations et le temps de calcul. Pour cela, un 
balayage de gauche a droite de 1 1 image est realise, sur 
chaque ligne de 1 ' image de reference ( c ' est -a-dire de 
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1' image obtenue avec eclairage en lumiere blanche) 
jusqu'au pixel qui a une intensite superieure a un 
seuil donne. Parmi les pixels obtenus pour toutes les 
lignes, celui avec la colonne d'ordre minimal a partir 
5 de la gauche determine la colonne extreme gauche, 
representee sch6matiquement sur la figure 7A par une 
ligne en traits interrompus 64. De la meme maniere, on 
determine la colonne extreme droite 66, la ligne 
extreme superieure 68 et la ligne extreme inferieure 

10 70. Le seuil applique (le meme pour tous les cotes de 
1' image) est defini selon le niveau de bruit de 
l f image. Ensuite, les pixels situes a l'interieur de 
l 1 image 62 sont codes a "O", tandis que les autres 
pixels sont codes a "1" . Ce codage est enregistre dans 

15 une matrice de meme taille que 1 1 image de reference. 

L'etape suivante permet de rechercher, pour 
chaque coeur, un point qui va representer le coeur tout 
entier, et auquel va etre affecte 1 ' intensite globale 

20 du coeur. Ce point peut etre notamment le barycentre de 
chaque coeur. 

La recherche est faite selon la methode de 
segmentation par croissance de region. Elle permet de 
regrouper les pixels contigus appartenant au meme 

25 coeur. Pour chaque pixel de 1 1 image de reference 62, on 
definit quatre connexions, selon le schema represents 
en figure 7C. Pour chaque connexion, un seuil T est 
defini (figure 10 : 6tape 80), tel que, pour tout pixel 
de coordonnee (xj, y^) : 

30 si I p (x 1 ,y 1 )>T, alors : p(x lf y 1 ) e 

ou Ip(xi, yi) designe 1' intensite du pixel de 
coordonnees (>^l/ Yl) et c k designe le coeur numero k. 
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A mesure que les pixels voisins sont 
rassembles, un numero k est attribue au coeur 
cor.respondant (figure 10 : £tape 82) . La methode de 
rassemblement s'applique d 1 une maniere recursive pour 
5 chaque pixel de 4 connexions. Le cas des pixels a moins 
de A connexions n'est constate que pour les pixels 
situes au bord de 1' image. Dans la recherche, il y a un 
test sur ces pixels et des qu'un pixel de bord est 
detecte, la recherche dans cette direction s'arrete et 
10 continue dans une autre direction. 

Une fois que les pixels d'un meme coeur sont 

rassembles, on limite (figure 10 : etape 84) le coeur 
par une fenetre, ou surface, dont : 

- le bord superieur est defini par le pixel p(x,y) doBt 
15 1 1 ordonnee y est maximale dans C^, 

- le bord inferieur est defini par le pixel p(x,y) dont 
1 * ordonnee y et minimale dans C^, 

- le bord lateral gauche est defini par le pixel p(x,y) 
dont I'abscisse x est minimale dans C^, 

20 - le bord lateral droit est defini par le pixel p(x,y) 
dont l ! abscisse x est maximale dans C^, 

II est ensuite possible de calculer le 
barycentre des pixels regroupes sous la designation 
"coeur n° k" (figure 10 : etape 86). On peut aussi 

25 determiner, parmi les pixels d'un coeur determine, le 
pixel d ! intensite maximale dans 1 ' image de reference. 

Le seuil T est determine a partir de 
1 1 histogramme cumule des intensites de 1 ' image de 
reference (en ne tenant compte que des pixels situes a 

30 l'interieur de la fibre, en utilisant le masque 
defini), et il est base sur les informations donnees 
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par le f abricant . De maniere plus precise, a partir des 
donn£es relatives au diametre du coeur, au nombre des 
coeurs de la fibre, le nombre moyen de pixels est 
calcule selon la relation suivante : 

N pc = (N p bxA mc ) /A mb 



ou A mc designe la surface d 1 un coeur (en pm 2 ) calculee 
a partir du diametre du coeur ; A mb designe la surface 

de la section de la fibre rnulticoeurs (en pm 2 ) definie 
selon les donnees du ; N pb designe le nombre de pixels 
dans 1' image de la fibre rnulticoeurs, ou encore la 
surface de la section de la fibre rnulticoeurs en pixels 
definie a partir de 1' histogramme cumule ; N pc designe 
le nombre de pixels du coeur, ou encore la surface 
moyenne de chaque coeur individuel, en pixels. 

II est possible que les fibres se deforment 
pendant leur etirage a chaud, et le diametre des coeurs 
peut done diminuer ou augmenter suivant une direction. 
Les coeurs n'ont alors plus la forme circulaire. Done, 
le nombre de pixels obtenus n'a -qu'une valeur moyenne. 
En multipliant le nombre de coeurs par la surface (en 
nombre de pixels) d ? un coeur individuel, on obtient la 
surface "utile" de la fibre rnulticoeurs. Ensuite, a 
partir de 1' histogramme cumule tel que celui illustre 
sur la figure 7D, et qui regroupe, en abscisse, le 
niveau cumule de gris et en ordonnee le nombre de 
pixels, on definit le niveau de gris qui correspond a 
cette surface "utile". L 1 histogramme donne, pour chaque 
classe de niveaux de gris, le nombre de pixels dont le 
niveau de gris se situe dans cette classe ou dans une 
classe de niveau de gris inferieur : le seuil de 
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segmentation T est defini par la valeur de niveau de 
gris qui correspond a la surface utile, mesuree en 
pixels. La recherche se fait en se deplacant de la 
droite vers la gauche (en niveau de gris), sur 
1' histogramme cumule, jusqu'a ce que la valeur de la 
surface utile (en nombre de pixels) des coeurs soit 
atteinte . 



10 



15 



II faut considerer de maniere particuliere le 
cas des coeurs fusionnes. En effet, au cours de la 
realisation des fibres multicoeurs, certains coeurs 
individuels peuvent se trouver suf f isamment proches 
l'un de 1' autre pour que, lors de l'etirage a chaud, 
ils fusionnent. L ■ image de ces coeurs individuels, sur 
1' image de reference ou sur une image objet, n'est 
alors plus de forme sensiblement circulaire, comme sur 
la figure 5. Au contraire, comme illustre sur la figure 
8, qui represente une partie d'une section d'une fibre 
multicoeur, les groupes de coeurs fusionnes 72, 74 
20 presenters une extension suivant une certaine 
direction. La reference 72 designe un groupe de deux 
coeurs fusionnes, la reference 74 un groupe de trois 
coeurs fusionnes. 

On peut realiser un test (figure 10 : etape 88) 
25 sur la surface de chaque coeur, obtenue par 
regroupement de la maniere deja expliquee ci-dessus, 
pour determiner si 1 ' on est en presence d'un groupement 
de coeurs fusionnes, ou pas. Ainsi, on peut definir un 
seuil sur le nombre de pixels de chaque classe C k 
obtenue pour determiner si on est dans le cas de la 
fusion de plusieurs coeurs (figure 10 : etape 92), ou 
non. Ce seuil est donne par la relation : 



30 
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N pf >(l+A n )N pc , 

ou N pf designe le nombre de pixels calcule dans un 
coeur fusionne (la signification de N pc a deja ete 

donnee ci-dessus) . Le parametre A n depend de 

5 1 1 importance de l'effet d ■ intermodulation (ou de 
couplage entre les coeurs, tel que deja defini plus 
haut) au niveau de l f image. En general, il varie entre 
0,3 et 0,8, et depend aussi de la valeur T de seuil 
globale. Par exemple pour 6000 coeurs, avec peu d'effet 

10 de couplage, une valeur elevee de A n (de 1 ' ordre 0,6- 

0,7) convient. Par contre, pour une fibre de 10000 
coeurs, ou l'effet de couplage est assez important, on 

choisit A n de 1' ordre 0,3-0,4. 

Une fois qu'un groupe de coeurs fusionnes est 
15 detecte, la question qui se pose est de savoir comment, 
sur une image obtenue a l'aide d'une fibre multicoeur, 
on peut retrouver la contribution de chacun des coeurs 
individuels qui sont fusionnes ensembles. A chaque 
groupe de coeurs fusionnes, tel que le groupe 74 
20 (figure 9A) on attribue une fenetre 76 qui englobe ce 
groupe. Cette fenetre est determinee, et passe, par les 
quatre points, ou pixels, dont une des deux coordonnees 
(colonne, ligne) est la plus grande ou la plus petite 
parmi tous les pixels du groupe 74. Ainsi, la fenetre 
25 76 passe par les points E et G qui sont respect ivement 
les points les plus a droite et les plus a gauche du 
groupement 74, et par les points F et H qui sont 
respectivement les points les plus hauts et les plus 
bas du groupement 74. 
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Le barycentre geometrique Pt>( x /Y) des pixels 

qui sont a l'interieur du groupement 74 est ensuite 
determine : on effectue la moyenne des nurr.eros de ligne 
et des numeros de colonne de :ous les pixels 
5 appartenant au groupement 74. 

On etablit ensuite la direction de la fusion 
(figure 10 : etape 94) . De maniere striate, cette 
direction (par exemple la direction A sur la figure 9A) 
peut en principe etre quelconque dans le plan de 
10 1' image. Prat iquement , on se fixe quatre directions 
dans la fenetre attribuee aux groupements fusionnes. 
Ainsi, sur la figure 9B, sont select ionnees les 

directions verticales D v , horizontale et les deux 

directions a ± 45° , D+ et D_ . Le nombre de pixels 

15 situes a l'interieur du groupement 74, le long de 
chacune de ces directions, est calcule. On retient, 
comme direction de fusion, celle de ces quatre 
directions le long de laquelle le nombre de pixels est 
maximum. Ainsi, sur la figure 9B, c'est la direction 

20 horizontale qui est retenue comme direction de 

fusion. La direction perpendiculaire (dans le cas de la 

figure 9B : D v ) est alors qualifiee de direction de 

coupure. 

On part ensuite du barycentre P^ de la classe 

25 74, et on effectue un deplacement, suivant la direction 
de fusion. Dans le cas ou il y a une fusion 
horizontale, le sens de deplacement est vers 1' extreme- 
gauche de la fenetre. Dans le cas d'une fusion 
verticale, le sens de deplacement est vers le haut. 

30 Dans le cas d'une fusion a +45° (-45° respecti vement ) 
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le sens de emplacement est vers 1' extreme droite, en 
haut (respectivement extreme gauche, en haut) . Ces sens 
respectent le balayage de 1 ' image (de haut en bas, et 
de gauche a droite) . On rencontre done successivement 
tous les pixels suivant la direction de fusion, et on 
mesure la variation d'intensite entre deux pixels 
adjacents. La figure 9C represente cette variation le 
long de la direction D h , de part et d' autre du 
barycentre P b . Lorsqu ' un minimum est rencontre, e'est- 
a-dire au point K, celui-ci est interprets comme un 
point de coupure de la fusion (figure 10 : etape 96). 
On fait done passer par le point K une droite de 
coupure, parallele a la direction de coupure D v . Cette 
droite de coupure delimite deux coeurs, 1 • un par 
rapport a 1' autre, a l'interieur d'un meme groupe de 
coeurs fusionnes. 

On peut ensuite, pour chacun des coeurs d'un 
meme groupe de coeurs fusionnes, calculer le 
barycentre, ou le point d'intensite maximale, auquel on 
attribue l'intensite totale transmise par le coeur. 



Lorsque la procedure de segmentation est 
terminee, on fait un test sur le nombre total des 

coeurs obtenus par segmentation. Si celui-ci est 

inferieur a un seuil donne (une fraction (par exemple 
80-85%) du nombre total de coeurs donne par le 
fabricant) 1' operation de segmentation est recommencee 
(figure 10 : etape 90), mais sans tenir compte des 
pixels qui appart iennent deja a un coeur detecte. Par 
ailleurs, au cours de cette deuxieme iteration, on 
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baisse le seuil de segmentation pour pouvoir regrouper 
les pixels de faible intensite. 

Le seuil defini pour cette operation depend du 
type de fibre et du nombre total de coeurs dans la 
fibre. En general, il est compris entre 70 et 80% du 
nombre total de coeurs donne par le fabricant. 

Ce cas se presente plutot pour les pixels qui 
se trouvent au bord de 1 1 image : a cause du bruit, ils 
sont de faible intensite. 



Une methode de determination d'un maillage 
triangulaire de 1 1 image va maintenant etre exposee. 
Cette methode permet d'attribuer, a chaque pixel, un 
triangle a partir duquel 1 * interpolation pourra etre 
15 realisee . 

Pour rassembler les coeurs de 1 ' image en groupe 
de trois, la methode la plus rapide est de determiner 
d'abord, pour chaque coeur, l 1 ensemble de ses voisins, 
puis de faire le regroupement a partir de cet ensemble. 
20 La matrice des points representatif s des coeurs 

est balayee par une fenetre carree, centree 'sur le plan 
representatif de chaque coeur. La taille de la fenetre 
depend de I'espace intercoeur et du type de fibre ; ces 
informations sont donnees par le fabricant de la fibre. 
Typiquement, des fenetres de taille 10 pixels x 10 
pixels, ou encore 7 pixels x 7 pixels, conviennent. 
Tous les coeurs dont le point representatif se situe a 
l'interieur de cette fenetre font partie du voisinage 
du coeur central de la fenetre carree. La figure 11A 

represente un point representatif d'un coeur, la 

fenetre 100 centree sur ce coeur et l 1 ensemble des 
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coeurs C k#1 ,...C kf5 qui font partie du voisinage du 
coeur C k . 

Les triangles sont ensuite construits a partir 
d'un coeur central et de deux coeurs voisins 
consecutifs. Ainsi, sur la figure 11A, le triangle 102 
est determine a partir des coeurs C k , C k)1 , C k/2 . Le 
triangle 104 est determine a partir du coeur Ck et de 
c k,2 et C k/3 . Une fois qu ' un triangle est determine, 
par exemple le triangle 102, tous les pixels qu'il 
contient sont affectes d'un codage, qui represente ce 
meme triangle. Autrement dit, pour chaque pixel a 
l'interieur du triangle dont les sommets sont C k , C k ,i, 

Ck,2/ 1 ' interpolation sera faite par rapport a ces 

trois sommets. 

De maniere plus precise, pour determiner les 

pixels appartenant a un triangle de sommet Si, S 2 , S 3 , 

une procedure est appliquee, qui va etre decrite en 
liaison avec la figure 11B : 

1. Determiner les sommets superieur Si, inferieur S 3 et 

le sommet "milieu" S 2 (situe entre les deux autres 
sommets, sur un axe vertical) . 

2. Diviser la surface du triangle en deux parties 
separees par la ligne horizontale Lh qui contient le 
sommet "milieu" S 2 du triangle. 

3. Determiner les cotes P i2 , P 13 du triangle qui lient 
le sommet superieur aux sommets inferieur S 3 et 
"milieu" S 2 . 
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4. Balayer ligne par ligne a partir du sonunet superieur 
Si et determiner la tranche de la ligne situee entre 

les deux cotes definis P12 et P13. 

5. Coder les pixels se situant sur la tranche de la 
5 ligne, et leur attribuer les trois coeurs Si, S2, S3 

constituant les sommets du triangle. 

6. A chaque traitement d'un pixel, verifier si il est 
code ou non (on evite ainsi le recouvrement des 
triangles) . 

10 7. Refaire les memes operations (3-6) pour la partie 
inferieure du triangle sauf que cette fois-ci on 

tient compte du cote P23 du triangle qui lie le 

sommet "milieu" S2 au sommet inferieur S3. 

De preference, avant de traiter un triangle, on 
15 determine si un des cotes du triangle est aligne 
horizontalement avec une ligne de pixels ou pas. Si le 
triangle n'est aligne avec aucune ligne de pixels, on 
procede de la maniere decrite ci-dessus (etapes 1 a 7). 
Si un des cotes est aligne avec une ligne de pixels, il 
20 n'est alors pas necessaire de faire la dif f erenciat ion, 
pour le triangle concerne, entre partie du haut et 
partie du bas du triangle. Le fait de faire ce test 
permet de diviser le temps de traitement par un facteur 
de I'ordre de 2. En effet, dans le.cas particulier d 1 un 
triangle dont un des cotes est aligne sur une ligne de 
pixels, le temps de remplissage du triangle avec les 
pixels est nettement inferieur au temps de remplissage 
d'un triangle dont les sommets sont disposes de maniere 
quelconque . 

La procedure decrite ci-dessus ne permet pas 
d'inclure tous les pixels dans un triangle. Une fois la 
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procedure de triangulation terminee, il est preferable 
de verifier s • il y a des pixels qui n'ont aucun 
triangle associe. Ceci est aise puisque, ainsi qu ' on 
l'a dit ci-dessus, des qu'un pixel est attribue a un 
triangle, on lui affecte un codage correspondant a ce 
triangle. Le cas de pixels n'ayant aucun triangle 
associe va etre decrit en liaison avec la figure lie. 
Sur cette figure, on a represents deux triangles C k , 
c kl, C k2 et C Lf C L1 , C L2 , determines de la maniere 
decrite ci-dessus. Un pixel p (x,y) est compris entre 
ces deux triangles, mais n'appartient a aucun. Une 
methode consiste alors a effectuer une recherche dans 
une direction specifique, par exemple la direction 
horizontale, pour trouver les pixels les plus proches 
appartenant a un triangle. Par exemple, dans le cas de 
la figure 11C, une exploration suivant la ligne passant 
par le pixel p(x,y) permet de determiner les pixels 
p(x-d 1/y ) et p(x+d 2 ,y) situes respectivement a une 
distance d x et a une distance d 2 a gauche et a droite 
de p(x,y). ces pixels appartenant chacun a un triangle, 
ils ont ete prealablement codes en fonction du triangle 
auquel ils appart iennent . Chaque pixel determine done 
un ensemble de trois sommets. Les distances du pixel 
p(x,y) aux six sommets sont calculees et ordonnees. On 
25 en deduit les trois coeurs, ou points representatif s 
des coeurs, les plus proches du pixel p(x,y). On 
associe ces trois coeurs a ce pixel : ainsi, sur la 
figure 11C, la distance de p(x,y) a chacun des sommets 
c k' C kl , c k2 , C L , C L1 , c L2 sera successi vement 
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calculee. On en deduira que p(x,y) est associe au 
triangle C k , C L1 , C k2 - 

De cette maniere, tous les pixels appart iennent 
a un triangle et un seul. Le fait de tenir compte des 
5 pixels n 1 appartenant a aucun triangle permet 
d'ameliorer 1 1 image , en particulier au bord, par 
rapport au cas ou on n'en tient pas compte. 



La reponse de chaque coeur au signal lumineux 

10 qu'il regoit est definie par la somme des intensites 
des pixels appartenant a ce coeur. Pour avoir une 
evaluation exacte de cette reponse, on associe a chaque 
coeur une surface, par exemple limitee par les 
mediatrices entre le coeur central (voir figure 7#) 

15 et chacun des coeurs voisins C^i, Cj^, •■■ c k5- 

maniere encore plus precise, on prend les mediatrices 
entre le barycentre (ou le centre, ou le point 
d'intensite maximale) et le barycentre (ou le centre, 
ou le point d'intensite maximale) de chacun des coeurs 

20 voisins. D'autres delimitations sont possibles, mais 
celle faite a l'aide des mediatrices permet de diviser 
la surface du faisceau de coeurs en polygones 
adjacents ; de plus, chaque pixel de 1 • image appartient 
ainsi a une, et seulement une, region C^. En outre 

25 cette methode de division est aisee a mettre en oeuvre 
pour les quelques milliers de coeurs dont on dispose. 

La delimitation- de ces coeurs fusionnes par 
rapport aux autres coeurs, individuels et non 
fusionnes, se fait de la meme maniere. 

30 

Les techniques de rehaussement de 1 ' image 
reconstitute, mises en oeuvre dans le cadre de la 
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presente invention, vont etre exposees ci-dessous. Le 
rehaussement du contraste est important lorsque la 
perception visuelle de 1 ' information dans 1 ■ image est 
limitee par de faibles variations de niveau de gris des 
differentes structures de 1' image. Le rehaussement 
modifie 1 ' apparence d • une image, de maniere a ce qu'un 
observateur, ou une machine, puisse extraire plus 
facilement une certaine information desiree. 

Les techniques de rehaussement decrites ci- 
dessous sont basees sur la redefinition de la valeur 
des pixels de 1 ' image originale, par modification de 
1 ' histogramme. Le rehaussement est obtenu en 
redeployant 1 * histogramme des intensites pour donner un 
poids egale a tous les niveaux de gris. Ainsi, on force 
des pixels de niveau d'intensite tres dense a occuper 
d'autres niveaux moins denses, et l'echelle des gris 
est elargie. Ceci a pour consequence une augmentation 
sensible du contraste. En fait, cela revient a 
attribuer a chaque pixel une nouvelle valeur de niveau 
20 de gris, par 1 ' intermediaire d'une fonction F de 
transformation : 



4 " Vp^y.q! q m> 



ou Ip represente la valeur transformer du pixel p, l p 
la valeur originale du pixel p, x,y les coordonnees du 
pixel p dans 1' image, et q lf ...q m des parametres de 
trans formation . 



Selon une premiere approche (dite globale) la 
fonction F est appliquee a chaque pixel, independamment 
de son voisinage, et elle est determinee a partir de 
1 ' his togramme de toute 1 ' image . 
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La section de 1 ' histogramme a "elargir" est 
determinee a partir d ' un seuil S fixe a priori. Les 
intensites dont le nombre de pixels depasse le seuil S, 
font partie de la section £ "elargir". La valeur du 
5 seuil S est estimee a l'oeil et depend du niveau du 
contraste de 1 1 image originale. 

Comme la figure 12 le represente, une nouvelle 
gamme dynamique est definie. Elle est limitee par les 

deux intensites extremes, I max et I m i n , de la section. 

10 Les L niveaux de gris de la gamme dynamique (par 
exemple : L=256) sont distribues sur la nouvelle garnme 

dynamique, afin de determiner le pas A q de 

quantification : 

A q=L/ dmax-Imin+1) 
15 Aux intensites appartenant a la gamme 

dynamique, on attribue une nouvelle valeur selon la 
relation suivante : 

Ip = dP - Imin> xA q 
Pour les intensites qui sont en dehors de la 
20 section d T histogramme a "elargir", la transformation 

est differente : Si l p < I min : l p = I min ' \- Si 

-£p > ^max : lp = ^max • 

Ces operations permettent de repartir les 
pixels des regions denses dans des regions moins denses 
25 et d'obtenir un histogramme elargi et uniforme, done un 
meilleur contraste. 

Selon une autre approche, dite locale, on 
prevoit des regions de voisinage centrees sur chaque 



JNSDOCID: <WO 9742600A1_I_> 



WO 97/42600 PCT7EP97/02053 



10 



15 



20 



32 

pixel, et la fonction F est calculee selon 
1 'histogramme cumule de la region consideree. 

La relation suivante decrit la fonction F : 

Ip=Lx (Hc Um (Ip) )/(N w xN w ) . 
ou L designe la gamme dynamique (e.g., 256), H cum (l p ) 
la valeur de 1 ' histogramme cumule a I'intensite I p , et 
"N w xN w la taille de la region (fenetre glissante) . 

Ainsi, la nouvelle valeur du pixel est obtenue 
en fonction des caracteristiques de son voisinage. Les 
resultats sont plus satisf aisants que dans le cas de 
l'approche globale, mais il y a quelques problemes : 

- rehaussement exagere du bruit, dans des regions 
homogenes ou 1 ' histogramme des intensites presente 
des pics, 

- presence d'artefacts, dans les frontieres entre les 
regions balayees par la fenetre glissante, 

temps de calcul assez important, du au traitement 
individuel de chaque pixel de 1' image. 



Une autre methode locale permet de resoudre le 
probleme du bruit. 

En general, dans les images medicales, il y a 
deux types de bruit : 

- le bruit du fond present dans la zone "noire" de la 
25 gamme dynamique de 1' image, 

- le bruit du aux fluctuations dans des regions 
homogenes de l'objet (l'organe). Ces variations 
d'intensite sont parfois interpreters comme la 
texture de l'organe. 

30 En fait, ces deux types de bruit se 

caracterisent par des variations de niveau de gris dans 
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une gamme etroite de 1 ' histogramme global de 1 ! image. 
Mais, dans certaines regions de 1' image, la 
contribution du bruit peut etre dominante. Apres avoir 
applique le rehaussement local de 1 ' histogramme dans 
5 ces regions, les fluctuations du bruit deviennent plus 
visibles. Ce probleme peut etre resolu, en limitant le 
niveau du rehaussement dans ces regions. 

Pour cela, on caracterise le rehaussement du 
contraste par la valeur de la pente de la fonction de 

10 transformation liant 1* intensite d ' entree a 1' intensite 
de sortie. Une pente de valeur 1 ne presente aucun 
rehaussement, et une valeur elevee de la pente 
correspond a un rehaussement important. Done, la 
limitation du rehaussement du bruit est obtenue en 

15 limitant la pente de la fonction de transformation. 

D'apres la relation ci-dessus, Ip est 
proportionnelle a H cum (histogramme cumule) . 

Etant donne que 1 1 histogramme des intensites 
n'est que la derivee de 1 ' histogramme cumule, la pente 
20 de la fonction de transformation a une intensite 
donnee, autrement dit le rehaussement du contraste, est 
proportionnelle a la valeur de 1 1 histogramme a cette 
intensite : 

d(I^)/d(I p ) = (L/(N w xN w ) )x(H orig (I p ) ) 
25 ou H or i g (I p ) est la valeur de 1 1 histogramme original a 
I'intensite I p . 

Done pour limiter la pente de la fonction, on 
limite 1 1 histogramme a une valeur constante (C^) . 

En fait cette valeur est un multiple du nombre 
30 moyen de pixels par niveau d'intensite. Ce multiple est 
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defini par une constante Cf, fixee a priori : 

c r=c f xL(N w xN w ) • 

La constante Cf est choisie selon les 
caracteristiques de 1' image, et elle est fixe pour 

toutes les regions. Par contre, la valeur C£ varie 

selon la taille de la fenetre. Par exemple, au bord de 
1' image, ou la taille de la fenetre est plus petite, la 

valeur calculee pour C£ est plus elevee que celle des 
autres regions . 

Apres avoir limite 1 1 histogramme, celui-ci est 
normalise, car la gamme d'intensites d'entree est egale 
a la gamme d'intensites de sortie. Done l f ensemble des 

pixels qui etaient au-dela de la valeur limite c^, est 

redistribue sur toute la gamme dynamique (Figure 12) . 

Si R est le nombre de pixels redistribues par 
intensite et S le nombre de pixels depassant la valeur 

C( : R=S/L. 

Ensuite, on redefinit 1 1 histogramme apres la 
normalisation, de la maniere su'ivante : 

- si H orig (I p ) > C £, H t (I p )=c^ 

- sinon : H t (Ip)=H or ig(Ip)+R / 

ou H t (I p ) est la valeur de 1 1 his togramme transforme 
pour 1 1 intensite I p . 

A partir de cet histogramme, on determine le 
nouvel histogramme cumule et la fonction de la 
transformation . 
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La solution aux problemes de la presence 
d'artefacts et de temps de calcul important consiste, 
dans une methode locale telle qu'exposee ci-dessus de 
5 maniere generale, a modifier le pas cu deplacement de 
la fenetre glissante (K.W. Leszczynski et S. Shalev, 
Image and Vision Computing, vol. 7, N° 3, p. 205-209, 
1989) . 

L* image est echant illonnee selon une trame 
10 definie (rectangulaire, hexagonale) . Les pixels de la 
trame constituent les centres de la fenetre glissante. 

La fonction F est calculee pour chaque fenetre. 
La valeur obtenue donne la nouvelle valeur 
d'intensite pour le pixel qui est le centre de la 
15 region. Pour les pixels qui ne font pas partie de' la 
trame d ' echantillonnage , il faut interpoler la valeur 
d'intensite des centres les plus proches . 

Cette interpolation peut etre faite dans deux 
directions. Mais pour diminuer l'erreur accumulee due. a 
20 1 1 interpolation bilineaire, on ne definit la trame 
d 1 echantillonnage que dans une direction (horizontale 
ou verticale) et la fonction d 1 interpolation est 
appliquee sur les pixels colineaires avec les pixels de 
la trame. 

25 Dans la relation suivante Ip/WQ-Ip.W! 

representent les valeurs transformees du pixel p(x,y), 
obtenues a partir des fonctions calculees dans les deux 
regions adjacentes centrees respect ivement aux pixels 

wq et vji . La valeur d'intensite interpolee lp nt du 

30 pixel p est calculee selon la relation suivante : 
T int T t _t t 

1 p =I p /W 0 +[ax(I p# w 0 - Ip^Wi)], ou d=y p -y w0 et a=d/N w . 
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Pour les pixels situes au bord de 1« image, il y 
a seulement une region de voisinage, done 
1' interpolation n'a pas eu lieu pour ces pixels. Le 
choix de la direction de 1 ' interpolation depend de la 
structure de 1' image. Si le contraste varie rapidement 
selon la direction verticale, 1 • interpolation 
horizontale introduit moins d'erreur que 

1 ' interpolation verticale . 

Cette methode est moins "lourde" en terme de 
nombres d' operations arithmetiques . Par exemple, pour 
la methode locale dans sa forme generale exposee ci- 
dessus, le nombre de multiplications et divisions est 
egal a environ 3x(NxN), tandis que, selon la version 
interpolee, ce nombre est 2x(NxN), (NxN) etant la 
15 taille de 1' image. 

Une application de I 1 invention a ete decrite 
dans le cadre medical. 

D'autres applications, dans d'autres domaines, 
notamment industriels sont egalement concernes, des 
lors qu'un objet est etudie et qu 1 une image de cet 
objet est realisee a 1 ' aide d'une fibre multicoeurs ou 
d'une multifibre. 
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REVE ND I CAT I ON S 

1. Procede de traitement d ' une image obtenue 
par. une fibre multicoeur ou une multifibre, comportant 
les etapes suivantes : 

5 l.pour chaque coeur de la fibre multicoeur ou de la 
multifibre, calcul d'une intensite, egale a la somme 
des intensites des pixels de 1 1 image qui sont situes 
a l'interieur d'une zone, predelimitee , de 
delimitation dudit coeur. 

10 2. calcul d'une intensite normalisee de chaque coeur, en 
divisant son intensite transmise, calculee comme dans 
l'etape 1, par 1' intensite transmise par ce meme 
coeur dans une image de reference. 
3. calcul d'une valeur interpolee de 1' intensite ^de 

15 chaque pixel de l f image, en fonction des intensites 

normalisees de trois coeurs d'un maillage 
triangulaire predetermine, et de la position relative 
du pixel par rapport a ces trois coeurs. 

2. Procede de traitement selon la revendica tion 
20 1, la zone predeterminee autour de chaque coeur etant 

definie sur 1 1 image de reference, cette derniere etant 
transposee sur 1 1 image objet. 

3. Procede de traitement selon 1 1 une des 
revendications 1 ou 2, la zone predelimitee autour de 

25 chaque coeur etant la surface delimitee par les 
mediatrices passant entre le barycentre dudit coeur et 
le barycentre de chacun des coeurs voisins. 

4. Procede de traitement selon 1 1 une des 
revendications 1 ou 2, la zone predelimitee autour de 

30 chaque coeur etant la surface delimitee par les 
mediatrices passant, dans 1 1 image de reference, entre 
le pixel d' intensite maximum de chaque coeur, et le 
pixel d' intensite maximum de chacun des coeurs voisins. 
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5. Procede selon l'une des revendications 1 a 

4, chaque pixel etant muni de 1 • identif ication des 
trois coeurs predetermines comme etant les plus 
proches . 

5 6. Procede selon l'une des revendications 1 £ 

5, comportant en outre une 'etape de rehaussement de 
1' image objet interpolee. 

7. Procede selon la revendication 6, le 
rehaussement etant un rehaussement global. 
10 8 - Procede selon la revendication 6, le 

rehaussement etant local, a I'interieur d'une fenetre 
glissante . 

9. Procede selon la revendication 8, la pente 
de la fonction de transformation etant limitee. 

10. Procede selon l'une des revendications 1 a 
9, comportant au prealable les etapes suivantes de 
traitement de 1 ■ image de reference : 

l r . recherche du barycentre, ou du pixel d'intensite 

maximale, de chaque coeur. 
2 l .calcul d'une zone de delimitation autour de chaque 
coeur . 

11. Procede selon la revendication 10 
comportant le calcul d 1 un masque numerique qui delimite 
la zone de 1 1 image a 1' inter ieur de laquelle la 
recherche du barycentre, ou du pixel d'intensite 
maximale, et le calcul de la zone de delimitation 
autour de chaque coeur doivent etre effectues. 

12 . Procede selon 1 1 une des revendications 10 
ou 11, la recherche du barycentre, ou de pixel 
d'intensite maximale etant faite apres avoir determine, 
pour chaque pixel de 1 ' image de reference, 1' ensemble 
des pixels qui l'entourent et qui sont d'intensite 
superieure a un certain seuil donne T. 
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13. Procede selon la revendication 12,. dans 
lequel est determine, pour chaque pixel de 1' image de 
reference, si il s'agit d'un pixel situe en bord 
d' image, ou pas. 

14. Proced£ selon la revendication 12, la 
determination des pixels d'intensite superieure a T 
evoluant lineairement , et changeant de direction des 
qu'un pixel d'intensite inferieure au seuil T est 
rencontre ou bien des qu 1 un pixel de bord est detecte. 

15. Procede selon 1 * une des revendication 10 a 

14, la zone de delimitation associee a chaque coeur 
etant obtenue en associant, au barycentre ou au pixel 
d'intensite maximale de chaque coeur, une surface 
limitee par les mediatrices entre ledit coeur (Cj^.et 
chacun des coeurs voisins (C^, C^..., C^s) 

16. Procede selon l'une des revendications 10 a 

15, 1 ■ etape 1" comportant en outre une sous-etape de 
traitement des coeurs fusionnes, comportant la 
determination de la presence de coeurs fusionnes et du 
barycentre, ou du pixel d'intensite maximale, de chacun 
des coeurs de chaque groupe de coeurs fusionnes. 

17. Procede selon la revendication 16, la sous- 
etape de traitement d'un groupe de coeurs fusionnes 
comportant : 

- la determination du barycentre geometrique des pixels 
faisant partie du groupe de coeurs fusionnes. 

- la determination d'une direction de fusion, d 1 une 
direction de coupure, et des points de coupure de 
groupes de coeurs fusionnes. 

- la coupure du groupe de coeurs fusionnes selon la 
direction de coupure, en chaque point de coupure. 

18. Procede selon la revendication 10, 
comportant en outre une etape prealable de 
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determination d ' un maillage triangulaire de 1 1 image de 
reference, tel que chaque pixel fasse partie d'un 
triangle. 

19. Procede selon la revendication 18, I'etape 
de determination d'un maillage triangulaire 
comportant : 

- une sous-etape de determination, pour chaque 
barycentre ou point d'intensite maximale, des 
barycentres ou points d'intensites maximales 
immediatement voisins, et de determination des pixels 
qui appartiennent a chaque triangle, forme par deux 
barycentres ou points d'intensite maximale, choisis 
voisins 1 1 un de 1 1 autre parmi les barycentres ou 
points d'intensite maximale immediatement voisins du 

15 barycentre ou point d'intensite maximale considere, 

et par ce dernier, 

- pour les pixels n'ayant aucun triangle attribue par 
la sous-etape precedente, une sous etape de : 

* recherche des deux pixels les plus proches 
appartenant a un triangle, 

* determination des trois barycentres, ou points 
d'intensite maximale, les plus proches du pixel 
considere, parmi les six barycentres ou points 
d'intensite maximale, constituant les somrnets des 
deux triangles auxquels appartiennent les deux 
pixels les plus proches du pixel considere. 

20. Procede de traitement d'une image de 
reference obtenue par une fibre multicoeur ou une 
multifibre, comportant les etapes suivantes : 
l f . recherche du barycentre, ou du pixel d'intensite 

maximale, de chaque coeur, 
2'. calcul d'une zone de delimitation autour de chaque 
coeur . 
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21. Procede selon la revendicat ion 20 
comportant au prealable le calcul d'un masque numerique 
qui delimite la zone de 1 1 image a 1 ' interieur de 
laquelle la recherche du barycentre, ou du pixel 

5 d'intensite maximale, et le calcul de la zone de 
delimitation autour de chaque coeur doivent etre 
ef f ectues . 

22. Procede selon I'une des revendica'tions 20 
ou 21, la recherche du barycentre, ou du pixel 

10 d'intensite maximale etant faite apres avoir determine, 
pour chaque pixel de 1 1 image de reference, 1 1 ensemble 
des pixels qui 1 1 entourent et qui sont d'intensite 
superieure a un certain seuil donne T. 

23. Procede selon la revendicat ion 22, dans 
15 lequel est determine, pour chaque pixel de 1 1 image de 

reference, si il s'agit d'un pixel situe en bo.rd 
d 1 image, ou pas . 

24. Procede selon la revendicat ion 22, la 
determination des pixels d'intensite superieure a. T 

20 evoluant lineairement , et changeant de direction des 
qu'un pixel d'intensite inferieure au seuil T est 
rencontre ou bien des qu'un pixel de bord est detecte. 

25. Procede selon l'une des revendicat ion 20 a 
24, la zone de delimitation associee a chaque coeur 

25 etant obtenue en associant, au barycentre ou au pixel 
d'intensite maximale de chaque coeur, une surface 
limitee par les mediatrices entre ledit coeur (C^) et 

chacun des coeurs voisins (C^, C^. . . , Cj^) 

26. Procede selon l'une des revendications 20 a 
30 25, l'etape l 1 comportant en outre une sous-etape de 

traitement des coeurs fusionnes, comportant la 
determination de la presence de coeurs fusionnes et du 
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barycentre, ou du pixel d'intensite maximale, de chacun 
des coeurs de chaque groupe de coeurs fusionnes. 

27. Procede selon la revendication 26, la sous- 
etape de traitement d'un groupe de coeurs fusionnes 

5 comportant : 

- la determination du barycentre geometrique des pixels 
faisant partie du groupe de coeurs fusionnes. 

- la determination d'une direction de fusion, d f une 
direction de coupure, et des points de coupure de 

10 groupes de coeurs fusionnes. 

- la coupure du groupe de coeurs fusionnes selon la 
direction de coupure, en chaque point de coupure. 

28. Procede selon la revendication 20, 
comportant en outre une etape prealable de 

15 determination d'un maillage triangulaire de 1 1 image de 
reference, tel que chaque pixel fasse partie d'un 
triangle . 

29. Procede selon la revendication 28, 1' etape 
de determination d'un maillage triangulaire 

20 comportant : 

- une sous-etape de determination, pour chaque 
barycentre ou point d'intensite maximale, des 
barycentres ou points d'intensites maximales 
immediatement voisins, et de determination des pixels 

25 qui appartiennent a chaque triangle, forme par deux 

barycentres ou points d'intensite maximale, choisis 
voisins 1 1 un de 1' autre parmi les barycentres ou 
points d'intensite maximale immediatement voisins du 
barycentre ou point d'intensite maximale considere, 

30 et par ce dernier, 

- pour les pixels n'ayant aucun triangle attribue par 
la sous-etape precedents une sous etape de : 
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* recherche des deux pixels les plus proches 
appartenant a un triangle, 

.* determination des trois barycentres, ou points 
d'intensite maximale, les plus proches du pixel 
5 considere, parmi les six barycentres ou points 

d'intensite maximale, constituant les sommets des 
deux triangles auxquels appartiennent les deux 
pixels les plus proches du pixel considere. 

30. Procede pour etablir une image medicale, 
10 comportant : 

- 1 ' etablissement d ' une image d'une partie d'un 
patient, a l'aide d'une fibre multicoeur, 

- le traitement de 1' image selon 1 ' une des 
revendications 1 a 19. 

31. Procede pour etablir un diagnostic medical, 
mettant en oeuvre un procede pour etablir une image, 
selon la revendicat ion 30, le diagnostic etant realise 
sur le lieu ou 1 ' image a ete etablie. 

32. Procede pour etablir un diagnostic medical, 
mettant en oeuvre un procede pour etablir une image, 
selon la revendicat ion 30, le diagnostic etant realise 
a distance du lieu ou 1 ' image a ete etablie. 

33. Dispositif pour traiter une image obtenue 
par une fibre multicoeur ou une multifibre, 

25 comportant : 

1- Des moyens pour calculer, pour chaque coeur de la 
fibre multicoeur ou de la multifibre, une intensite 
transmise egale a la somme des intensites des pixels 
de 1' image qui sont situes a 1 ' interieur d'une zone, 
predelimitee, de delimitation dudit coeur. 
2. Des moyens pour calculer une intensite normalisee de 
chaque coeur, a partir de son intensite transmise, 
calculee comme dans 1 ' etape 2, par 1' intensite 



20 



30 



3NSDOCID: <WO 9742600A 1_l_> 



WO 97/42600 



44 



PCT/EP97/02053 



transmise par ce meme coeur dans une image de 
reference . 

3. Des moyens pour le calcul d'une valeur interpolee de 
1' intensity de chaque pixel de l f image, en fonction 
des intensites normalisees de trois coeurs d r un 
maillage triangulaire predetermine, et de la position 
relative du pixel par rapport a ces trois coeurs. 

34. Dispositif selon la revendication 33 
comportant en outre des moyens pour transposer 1' image 
de reference sur chaque image objet. 

35. Dispositif selon 1 1 une des revendications 
33 ou 34, comportant en outre des moyens pour munir 
chaque pixel de 1 1 identification des trois coeurs 
predetermines comme 6tant les plus proches . 

36. Dispositif selon l'une des revendications 
33 a 35, comportant en outre des moyens pour realiser 
un rehaussement de 1' image objet interpolee. 

37. Dispositif selon l'une des revendications 
33 a 36, comportant : 

1*. Des moyens pour rechercher un barycentre, ou un 
pixel d'intensite maximale, de chaque coeur d'une 
image de reference. 

2 1 . Des moyens de calcul d'une zone de delimitation 
autour de chaque coeur. 

38. Dispositif selon la revendication 37, 
comportant des moyens de calcul d f un masque numerique 
qui delimite une zone d'une image. 

39. Dispositif selon l'une des revendications 
37 ou 38, comportant des moyens pour determiner, pour 
chaque pixel de 1 1 image de reference, 1' ensemble des 
pixels qui I'entourent, et qui sont d'intensite 
superieure a un certain seuil donne T. 
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40. Dispositif selon la revendicat ion 39, 
comportant des moyens pour determiner, pour chaque 
pixel d'une image de reference, le nombre de pixels 
immediatement voisins. 
5 41. Dispositif selon 1 1 une des revendicat ions 

37 a 40, comportant en outre des moyens pour determiner 
la presence de coeurs fusionnes, et determiner le 
barycentre, ou pixel d'intensite maximale, de chacun 
des coeurs de chaque groupe de coeurs fusionnes. 
10 42. Dispositif selon la revendicat ion 41, 

comportant : 

- des moyens pour determiner le barycentre geometrique 
des pixels faisant partie d'un groupe de coeurs 
fusionnes. 

15 - des moyens pour determiner une direction de fusion, 
une direction de coupure, et des points de cpupure-.de 
groupes de coeurs fusionnes. 

- des moyens pour determiner la coupure du groupe de 
coeurs fusionnes selon la direction de coupure, en 

20 chaque point de coupure. 

43, dispositif selon la revendicat ion 37, 
comportant en outre des moyens de determination d'un 
maillage triangulaire d ! une image de reference. 

44. Dispositif pour 1 1 etablissement d'une image 
25 medicale, comportant une fibre multicoeur et un 

dispositif selon 1 1 une des revendications 33 a 43.. 
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